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Abstract 

The Relativistic Hartree Bogoliubov (RHB) model is applied in 
the description of ground-state properties of proton-rich odd-Z nuclei 
in the region 53 < Z < 69. The NL3 effective interaction is used 
in the mean-field Lagrangian, and pairing correlations are described 
by the pairing part of the finite range Gogny interaction D1S. The 
model predicts the location of the proton drip-line, the ground-state 
quadrupole deformations and one-proton separation energies at and 
beyond the drip-line, the deformed single-particle orbitals occupied by 
the odd valence proton, and the corresponding spectroscopic factors. 
The results of fully self-consistent RHB calculations are compared with 
available experimental data, and with predictions of the macroscopic- 
microscopic mass model. 
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1 Introduction 



Models based on the relativistic mean-field approximation provide a micro- 
scopically consistent, and yet simple and economical description of the nu- 
clear many-body problem. By adjusting just a few model parameters: cou- 
pling constants and effective masses, to global properties of simple, spherical 
and stable nuclei, it has been possible to describe in detail a variety of nu- 
clear structure phenomena over the whole periodic table, from light nuclei 
to superheavy elements When also pairing correlations are included in 
the self-consistent Hartree-Bogoliubov framework, the relativistic mean-field 
theory can be applied to the physics of exotic nuclei at the drip-lines. 

The relativistic Hartree-Bogoliubov (RHB) model represents a relativistic 
extension of the Hartree-Fock-Bogoliubov (HFB) framework. A unified and 
self-consistent description of mean-field and pairing correlations is especially 
important in applications to drip-line nuclei. At the drip-lines the separation 
energy of the last nucleons becomes extremely small or negative, the Fermi 
level is found close to particle continuum, and the lowest particle-hole or 
particle-particle modes are embedded in the continuum. The RHB model 
with finite range pairing interactions has been applied in studies of the halo 
phenomenon in light nuclei 0, [3] , properties of light nuclei near the neutron- 
drip [|]], reduction of the spin-orbit potential in nuclei with extreme isospin 
values 0, ground-state properties of Ni and Sn isotopes ||, and the defor- 
mation and shape coexistence phenomena that result from the suppression 
of the spherical N=28 shell gap in neutron-rich nuclei J7|. In particular, in 
Ref. || we have applied the RHB model to calculate properties of proton- 
rich spherical even-even nuclei with 14<Z<28 and N=18,20,22. It has been 
shown that RHB provides a very good description of binding energies, two- 
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proton separation energies and proton rms radii. Model predictions for the 
proton drip-line have been found in excellent agreement with shell-model 
calculations and with results of non-relativistic HF and HFB studies. The 
position of the two-proton drip-line for even-even nuclei with 10<Z<82 has 
been recently determined with the RMF+BCS model ||. Results computed 
with the NL3 effective interaction ]IUj have been compared with those of the 
Hartree-Fock+BCS model (effective force Skyrme SIII) [11|, the finite-range 
droplet model (FRDM) Jl2|, [R|, and with the available experimental infor- 



mation. In a very recent work [jrjj we have used the RHB model with finite 
range pairing for an analysis of the structure of deformed proton-rich odd- 
Z nuclei in the region 59<Z<69. The location of the proton drip-line, the 
properties of proton emitters beyond the drip-line, and the deformed single- 
particle orbitals occupied by the odd valence proton have been compared 
with available experimental data. 

The structure of proton-rich nuclei displays many interesting properties 
which are important both for nuclear physics and astrophysics. These nu- 
clei are characterized by exotic decay modes, such as the direct emission of 
charged particles from the ground-state, and /3-decay with large Q-values. 
The properties of most proton-rich nuclei should also play an important role 
in the process of nucleosynthesis by rapid-proton capture. The decay by 
direct proton emission provides the opportunity to study the structure of 
systems beyond the drip-line. Proton radioactivity from the ground-state 
is determined by the Coulomb and centrifugal terms of the effective poten- 
tial. For Z<50, nuclei beyond the proton drip-line exist only as short lived 
resonances, and ground-state proton decay probably cannot be observed di- 
rectly PI. On the other hand, the relatively high potential energy barrier 
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enables the observation of ground-state proton emission from medium-heavy 
and heavy nuclei. At the drip-lines proton emission competes with f3 + decay; 
for heavy nuclei also fission or a decay can be favored. The proton-drip line 
has been fully mapped up to Z=21, and possibly for odd-Z nuclei up to In 
[ 15 ] . Detailed studies of ground-state proton radioactivity have been reported 
for odd-Z nuclei in the two spherical regions 51<Z<55 and 69<Z< 83. The 
systematics of spectroscopic factors is consistent with half-lives calculated in 



the spherical WKB or DWBA approximations [16, 17, 18]. Recently reported 



proton decay rates [19| indicate that the missing region of light rare-earth 
nuclei contains strongly deformed systems at the drip-lines. The lifetimes of 
deformed proton emitters should provide direct information on the Nilsson 
configuration occupied by the odd proton, and therefore on the shape of the 
nucleus. 

In the present work we extended the analysis of Ref. [14], and report a 
detailed study of deformed odd-Z proton-rich nuclei with 53<Z<69. In Sec. 
2 we present an outline of the relativistic Hartree-Bogoliubov theory and 
discuss the various approximations and symmetry constraints. In Sec. 3 the 
model is applied in the calculation of ground-state properties of proton-rich 
isotopes at and beyond the drip-line. The results are summarized in Sec. 4. 



2 Relativistic Hartree Bogoliubov model with 
finite range pairing interaction 

2.1 The Relativistic Mean Field Model 

The relativistic mean field theory presents a phenomenological description 
of the nuclear system. The theory is based on very simple concepts: (i) 
nucleons are described as point particles, (ii) these particles strictly obey the 
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rules of relativity and causality, and (iii) they move independently in mean 
fields which originate from the nucleon-nucleon interaction. Conditions of 
causality and Lorentz invariance impose that the interaction is mediated by 
the exchange of point-like effective mesons, which couple to the nucleons 
at local vertices. The model parameters are the meson masses and their 
coupling constants. The values of the parameters are adjusted to properties 
of nuclear matter and experimental data on finite nuclear systems. 

On this phenomenological level the different mesons are classified accord- 
ing to the quantum numbers of spin, parity and isospin. The concepts of 
simplicity, and the desire to keep the number of phenomenological parame- 
ters as small as possible, impose that only the most important meson fields 
should be included in the model. The experimental fact that parity is con- 
served to a very high degree of accuracy, on the mean field level automatically 
excludes mesons with negative intrinsic parity. Since the model describes a 
system with two types of particles: neutrons and protons, it is necessary to 
introduce both isoscalar and isovector meson fields. Most versions of the 
relativistic mean-field model, however, do not contain the scalar isovector 
meson 5. From the presently available experimental data it is difficult to 
determine its effective coupling constant and effective mass, i.e. the same 
level accuracy is obtained with or without the inclusion of the 5-meson. On 
the other hand, essential for a correct description of properties of finite nuclei 
is the vector isovector p-meson. 

The nucleons are described as relativistic particles moving independently 
in local mean-field meson potentials. The single-nucleon relativistic dynam- 
ics is described by the Dirac equation. In contrast to the Schrodinger equa- 
tion, which contains a single mean-field potential of Woods-Saxon shape and 
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an independently adjusted spin-orbit potential, the Lorentz structure of the 
Dirac equation distinguishes between the scalar potential S(r), which results 
from the exchange of spin-zero mesons, and a vector potential induced by 
spin-one mesons. In addition, the spin-orbit potential is automatically in- 
cluded in the equation of motion, i.e. it does not require an independent 
strength parameter. The vector potential (Vo(r), V(r)) is a four-vector un- 
der Lorentz transformations. It is similar in structure to the electromagnetic 
potentials (v4 (r), A(r)) which occur when the Dirac equation is applied to 
systems in atomic physics. The time- like component of the vector field Vo(r) 
corresponds to the Coulomb field A (r), and the three space-like components 
V(r) are equivalent to the magnetic potential A(r). 

In calculations of properties of ground-state proton emitters we neglect 
the space-like components of the vector fields, i.e. we neglect nuclear mag- 
netism. For nuclei with even number of protons and neutrons this is certainly 
justified by the time-reversal invariance of the mean field. In nuclei with an 
odd number of protons or neutrons, time-reversal symmetry is broken. The 
odd particle induces polarization currents and time-odd components in the 
meson fields. The time-odd components are essential for a description of 
magnetic moments [ 2(| , and of moments of inertia in rotating nuclei f2l| , ^2 



The effect on deformations and binding energies, however, is very small and 
can be neglected to a good approximation. We are then left with the time- 
like component V(r) (for simplicity we omit the subscript henceforth). The 
Lorentz structure of the theory implies that the scalar field is attractive, and 
that the time-like component of the vector fields is repulsive. Hence, the 
scalar field can be used to characterize in an economic way the attractive 
part of the nucleon-nucleon interaction at intermediate distances, and the 
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vector fields describe the repulsion at short distances. 

Under these simplifying assumptions the single-nucleon dynamics is de- 
scribed by the Dirac equation, which contains only the time-like component 
of the vector potential V and the scalar potential S. In the stationary case 
it reduces to an eigenvalue problem 

h D iPi = {ap + V(T) + P[m-S(T)]}iPi= enPi. (1) 

This equation contains the four- dimensional Dirac matrices a and f3, and the 
nucleons with rest mass m are represented by the four-dimensional vectors 
ipi. Expressed in terms of the Pauli spin matrices, the equation reads 

( m %> +v - m + p s+v){ f 9 )r £t { f 9)- (2) 

The single-particle wave functions ipi are four- dimensional spinors. The sub- 
script i denotes the quantum numbers that specify the single-particle state, 
and is the corresponding energy. In the no-sea approximation the Dirac 
sea of solutions with negative eigenvalues does not contribute to the densities 
and currents. Performing only the summation over the occupied orbitals in 
the Fermi sea, the bilinear products of wave functions are used to construct 
two types of densities: (i) the baryon density 

A A 

p{v) = £>M = Eltf to/ito +9H*)9M, (3) 

i=l i=l 

which is the zero component of the four- vector baryon current, and (ii) the 
scalar density 

A A 

P.(r) = E^ = TtftWiW -9H*)gM- (4) 

i=i i=i 

The fields V(r) and S(r) are computed by averaging over the interactions 
induced by the exchange of vector and scalar mesons, respectively 

V(r) = J v v (r,r')p(r')d 3 r, (5) 



S(r) = |^(r,r')p s (r')rf 3 r. (6) 

The two-body interactions v s (r, r') and v v (r, r') are of Yukawa type. They 
correspond to the exchange of the scalar isoscalar a-meson, and vector mesons: 
uo (isoscalar) and p (isovector). The resulting fields are defined 

S(r) = -g a v(*), (7) 
V(r) = 9u u(t) + g pP3 (T) + Ao(^ ( 8 ) 

where cr(r), u;(r), £3(1") are the classical mesons fields, and A (r) is the 
Coulomb field (exchange of photons). The mesons fields are solutions of 
the Klein-Gordon equations 

(-A + m a )a(r) = -g aPs (r), (9) 

(-A + m u )u(r) = g u p(r), (10) 

(-A + m p )p 3 (r) = g p (p n (r) -p p (r)), (11) 

-AA (r) = e 2 p c (r), (12) 

m CT , and m p are the masses of the a-meson, the a;-meson, and the p- 
meson, respectively. g a , g^, and g p are the corresponding coupling constants 
for the mesons to the nucleon, and e 2 /47r = 1/137.036. p c (r) is the charge 
density. The interactions that result from the Yukawa and Coulomb type 
Green's function are 

„2 -m a \r-r'\ 

v s (r,r') = -§z--. 7r , (13) 

47r |r — r | 

2 -rriwlr-r'l fl 2 _ m | r _ r '| 2 i 

wV ' ' 4tt |r-r'| 4tt |r-r'| 4?r |r - r'| ' V ; 

where r are the isospin matrices. 

In the mean-field approximation the mesons are introduced only on a phe- 
nomenological level as carriers of the fields. On a more fundamental level, 



one should certainly try to relate these effective mesons to physical mesons. 
For the vector mesons u and p this relation is relatively well established. The 
effective masses computed from the experimental properties of finite nuclei 
are very close to the free values. In fact, in many successful parameter sets 
reported in the literature, one finds that the bare values have been used for 
the masses of the u and p mesons. This choice reduces the number of free 
model parameters by 2. In the scalar sector, however, the phenomenological 
level of the model becomes evident. A physical scalar isoscalar meson has 
not been observed. On the other hand, it is well known that the one-pion ex- 
change determines the nucleon-nucleon interaction at large distances. In the 
region of intermediate distances the interaction is predominantly the corre- 
lated two-pion exchange and the two-pion exchange with a A-particle in the 
intermediate state. Both processes result in a parity-conserving mean field, 
which in the RMF model is described by the scalar field S(r). In princi- 
ple, the Lorentz structure would allow a pseudo-scalar field which originates 
from the one-pion exchange. On the Hartree level, however, this contribution 
vanishes due to parity conservation, and therefore the pion contributes only 
through the two-pion exchange. 

Summarizing in a formal way our considerations, RMF can be formulated 
as an effective field theory, starting from the Lagrangian density originally 
introduced by Walecka |23| : 



L = if) (z*7 • d — m) ip + -(da) 2 — -m a a 2 



(1 — r 3) 

- g a ipaip - g^ipj ■ utp - g p tpj ■ pr^ - e^7 ■ A ip .(15) 

Vectors in isospin space are denoted by arrows, and bold-faced symbols will 
indicate vectors in ordinary three-dimensional space. Q^, R^, and are 
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the field tensors of the vector fields u, p, and of the photon: 

QT = (16) 
RT = d»p u -d u p» (17) 
F n» = d tx A v -d v A li . (18) 

If the bare masses m, m w , and m p are used for the nucleons and the u and 
p mesons, there are only four free model parameters: m a , g a , g^ and g p . 
Their values can be adjusted to the experimental data on just few spherical 
nuclei. This simple model, however, is not flexible enough for a quantitative 
description of properties of complex nuclear systems. An effective density 
dependence is introduced [|4]] by replacing the quadratic cx-potential \m 2 r a 2 
with a quartic potential U (a) 

U(a) = l m la 2 + ±g 2 a 3 + ^g 3 a\ (19) 

The potential includes the nonlinear a self-interaction, with two additional 
parameters g 2 and #3. The corresponding Klein-Gordon equation (|9]) becomes 
nonlinear, with a cr-dependent mass m 2 (a) = m 2 + g 2 a + g^u 2 . More details 
on the relativistic mean- field formalism can be found in Refs. |T|, [531, |25], pEfl . 

The essential structure of the model is based on relativity. In a descrip- 
tion of the ground-state properties of nuclei we are certainly faced with the 
question: why is it necessary to use a relativistic formulation? In fact, the 
kinetic energies and the Fermi momenta are relatively small compared to 
the rest mass of the nucleon. Therefore, relativistic kinematics can be ne- 
glected. The Dirac equation, however, contains more than just the difference 
V(r) — S(r). Compared to an equivalent Schrodinger equation with a poten- 
tial well of ~50 MeV (small compared to the nucleon rest mass of 939 MeV), 
the Dirac equation contains two potentials: V(r) and —S(r), which are both 
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large (~ 350 and ~ —400 MeV) and opposite in sign. Therefore, one needs 
a relativistic framework in order to properly describe the dynamics of these 
strong potentials. In the equation (||) for the large component /, only the 
difference between of V and S enters, and this potential is, in fact, small 
compared to the nucleon rest mass. In the lower equation for the small com- 
ponents g, the sum V+S occurs. This term is large and it produces the strong 
spin-orbit potential, which is essential in nuclear structure calculations. One 
of the advantages of working in a relativistic framework is also being able to 
determine the strength and the shape of the spin-orbit term self-consistently. 
This is especially important for a correct description of spin-orbit splittings 
in regions of nuclei far from the valley of (3 stability, where the extrapolation 
of effective strength parameters becomes questionable. Another example of 
the importance of describing the single-nucleon wave functions with Dirac 
spinors, is the occurrence of approximate pseudo-spin symmetry in nuclear 
spectra. The symmetry has its origin in the approximate cancellation of the 



scalar and vector potentials in the single-nucleon Dirac equation [2y| . 

It should also be emphasized that the smallness of V — S leads to relatively 
small Fermi momenta and allows, in principle, a non-relativistic reduction of 
the Dirac equation to a Schrodinger equation with momentum-dependent 
potentials. The resulting non-relativistic theory with additional spin- and 
momentum-dependent terms requires, in general, more adjustable parame- 
ters, and one can argue that its predictive power is reduced as compared to 
a fully relativistic theory. 

2.2 Relativistic Theory with Pairing Correlations 

The inclusion of pairing correlations is essential for a correct description of 
structure phenomena in spherical open-shell nuclei and in deformed nuclei. 
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The relativistic mean field theory, as described in the previous section, does 
not include such correlations. The Dirac Hamiltonian contains only single- 
particle field operators with the structure i/j^i/j; pairing correlations can only 
be described in a generalized single-particle theory by field operators ip^ip^ 
and if)if), which do not conserve the particle number. Pairing correlations 
are therefore often included in a phenomenological way with the simple BCS 



approximation ||28|| . However, the BCS model presents only a poor approx- 
imation for exotic nuclei far from the valley of /3-stability. The physics of 
drip-line nuclei necessitates a unified and self-consistent treatment of mean- 
field and pairing correlations. In the non-relativistic approach, properties of 
drip line nuclei are consistently described in the framework of Hartree-Fock- 
Bogoliubov (HFB) theory [BjJ. The relativistic extension of the HFB theory 



was introduced in Ref. [30]. The starting point is again the Lagrangian 
(p"5*|). By performing the quantization of the meson fields, the relativistic 
Hartree-Bogoliubov equations are derived using Green's function techniques. 

As in non-relativistic HFB theory, the ground state of a nucleus |$ > is 
described as the vacuum with respect to independent quasi-particle opera- 
tors, which are defined by a unitary Bogoliubov transformation of the single- 
nucleon creation and annihilation operators. The generalized single-nucleon 
Hamiltonian contains two average potentials: the self-consistent Hartree- 
Fock field T which encloses all the long range ph correlations, and a pairing 
field A which sums up the ^^correlations. In the relativistic extension the 
Hartree approximation is employed for the self-consistent mean-field, and the 
resulting Relativistic Hartree-Bogoliubov (RHB) equations read 

(ho-m-X A \(U k (r)\ (U u (r)\ 

{ -A* -h D + m + X ) { VUr) )- Eh [ V h (x) ) (20) 

where hp is the single-nucleon Dirac Hamiltonian (p]), and m is the nucleon 
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mass. The chemical potential A has to be determined by the particle num- 
ber subsidiary condition, in order that the expectation value of the particle 
number operator in the ground state equals the number of nucleons. The 
column vectors denote the quasi-particle wave functions, and E k are the 
quasi-particle energies. 

As it is well known, from the non-relativistic Hartree-Fock-Bogoliubov 
theory and the detailed discussion in Ch. 7 of Ref. f3lH , for each of the 
positive quasiparticle energies E k > with eigenvector (U k ,V k ) exists an 
eigensolution — E k with eigenvector (V k *,U k ). One has to choose for each 
quantum number k one of these solutions. In Hartree-Fock calculations this 
corresponds to the choice of a certain occupation-pattern for the Slater de- 
terminant. For the calculation of the ground state of an even-even nucleus 
one occupies the lowest levels starting from the bottom of the well, i.e. one 
chooses the lowest A single-particle levels. In the HFB case the ground- 
state of an even-even system is obtained by choosing only the positive quasi- 
particle energies (see Ref. |JT| for exceptions of this rule in the case of rotating 
nuclei). To describe odd nuclei one has to block at least one level ho. This 
means that one has to choose for ko the negative eigenvalue — E ko with the 
eigenvector (V^, U ko ). Therefore, in the following equations summations over 
k run only over half of the eigensolutions of the RHB-equation ( p0|) chosen 
according to this prescription. 

The RHB equations are solved self-consistently, with potentials deter- 
mined in the mean-field approximation from solutions of Klein-Gordon equa- 
tions 



A + m 2 a a{r) 



9aps(r) -g 2 (J 2 (v) -g 3 a 3 (r) 



(21) 



(22) 
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-A + mJJ p°(r) = g p p 3 (r) (23) 
—AA°(r) = e Pp (r), (24) 

for the sigma meson, omega meson, rho meson and photon field, respectively. 
Due to charge conservation, only the 3rd-component of the isovector rho 
meson contributes. The source terms in equations (|2l"D to (|24"D are sums of 
bilinear products of baryon amplitudes 

Ps(r) = EnWm (25) 

k 

Pv(r) = ^V k \r)V k (r), (26) 

k 

Pa(r) = J2Vi(r)r 3 V k (r), (27) 

k 

Pe m (r) = E#)^ (r) ' (28) 

where the sums run over all the all the quasiparticle states with positive 
energy in even-even nuclei. In odd systems at least one level ko has to be 
blocked, i.e. the corresponding eigenvector V k (r) has to be replaced by 

The pairing field A in ( pUP is an integral operator with the kernel 

A a6 (r, r') = ^ E Vabcdir, r')K cd {r, r'), (29) 

c,d 

where a, b, c, d denote quantum numbers which specify the Dirac indices of 
the spinors, V a b c d(r, r') are matrix elements of a general two-body interaction. 

The eigensolutions of Eq. fl2"U|) form a set of orthonormal single quasi- 
particle states. The corresponding eigenvalues are the single quasi-particle 
energies. The self-consistent iteration procedure is performed in the basis of 
quasi-particle states. As shown by Bloch and Messiah in |3"2"fl , any Hartree- 
Fock-Bogoliubov wave function can also be represented in the form of a BCS 
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function. For systems with an even number of particles 

m = n(«*+v*M)i->> ( 3 °) 

and for odd particle number 

i $ mo) = < n k+vW)i->> ( 3i ) 

where /j,q is the level occupied by the odd particle (blocked level). 

The operators aj, and ajj create particles in the canonical basis, and the 
occupation probabilities are given by 

" 2 V - m " A ) 2 + A ' 

£ M = ()u|/i£>|ju) and A^ = um|A|/Z^ are the diagonal elements of the Dirac 
single-particle Hamiltonian and the pairing field in the canonical basis, re- 
spectively. In contrast to the BCS framework, however, neither of this fields is 
diagonal in the canonical basis. The basis itself is specified by the requirement 
that it diagonalizes the single particle density matrix p(r, r') = ^4( r )^fc ( r 0- 

k 

The transformation to the canonical basis determines the energies and oc- 
cupation probabilities of single-particle states, which correspond to the self- 
consistent solution for the ground state of a nucleus. 

According to the Bloch- Messiah theorem the canonical basis is connected 
to the quasiparticle basis through a unitary transformation (usually called 
C-transformation). Only in the HF+BCS case this transformation is pro- 
portional to unity. Therefore, the blocked level /xq i n Eq. Q5TD, which is a 
basis vector of the canonical basis, is in general different from the blocked 
quasiparticle level k discussed after Eq. fl20D . However, according to the 
Bloch- Messiah theorem a density matrix p(r, r'), obtained after blocking the 
quasiparticle level ko, has always one eigenvalue 1 and one eigenvalue in 
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addition to the doubly degenerate eigenvalues v 2 . Neglecting degeneracies, 
there is a unique relation between the quasiparticle level ko and the blocked 
level /io in the canonical basis. In case of symmetries both levels /io and ko 
have identical quantum numbers while in cases where the matrix C is close 
to unity, the corresponding wavefunctions are similar. 

In many applications of the relativistic mean-field model pairing correla- 
tions have been taken into account in a very phenomenological way in the 
BCS model with the monopole pairing force, adjusted to the experimental 
odd-even mass differences. This framework obviously cannot be applied to 
the description of the coupling to the particle continuum in nuclei close to 
drip-line. The question therefore arises, which pairing interaction V a b c d{r, r') 



should be used in Eq. (18). In Ref. |30| a fully relativistic derivation of the 
pairing force has been developed, starting from the Lagrangian (1). Using 
the Gorkov factorization technique, it has been possible to demonstrate that 
the pairing interaction results from the one-meson exchange (a-, u- and p- 
mesons), and is therefore of the form fllTf ) and In practice, however, it 
turns out that the pairing correlations calculated in this way, with coupling 
constants taken from the standard parameter sets of the RMF model, are too 
strong. The repulsion produced by the exchange of vector mesons at short 
distances results in a pairing gap at the Fermi surface that is by a factor 
three too large. However, as has been argued in many applications of the 
Hartree-Fock-Bogoliubov theory, there is no real reason to use the same ef- 
fective forces in both the particle-hole and particle-particle channels. In the 
ladder approximation, the effective interaction contained in the mean-field f 
is a G-matrix, the sum over all ladder diagrams. The effective force in the 
pp channel, i.e. in the pairing potential A, should be the K matrix, the sum 
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of all diagrams irreducible in the pp-direction (see for instance |P5fl). Since 
very little is known about this matrix in the relativistic approach, in most 
applications of the RHB model a phenomenological pairing interaction has 
been used, the pairing part of the Gogny force, 

V pp (l,2) = e-« ri - r2)UH)2 (Wi + B i P' 7 -HiP T -M^P^, (33) 

i=l,2 



with the set D1S |34| for the parameters fa, Wi, B i} Hi and Mj (i = 1, 2). 

Although this procedure formally breaks the Lorentz structure of the 
equations, one has to keep in mind that pairing in itself is a completely non- 
relativistic phenomenon. The essential point in a relativistic description of 
nuclei is the difference between the scalar and the vector density, which yields 
saturation, together with the approximate cancellation of the scalar attrac- 
tion by the vector repulsion, for the large components in the Dirac equation. 
The addition of the scalar and vector potentials for the small components, re- 
sults in large energy splitting between spin-orbit partners. The admixture of 
small components through the kinetic term up in the off-diagonal part of the 
Dirac equation, has no counterpart in the pairing field. Relativistic effects 
are only important for the mean field part of Hartree- Bogoliubov theory. 
The pairing density k and the corresponding pairing field, result from the 
scattering of pairs in the vicinity of the Fermi surface. The pairing density 
is concentrated in an energy window of less than 50 MeV around the Fermi 
level, i.e. the contributions from the small components to the pairing tensor 
K are extremely small. The word relativistic in RHB, therefore, applies al- 
ways the the Hartree p/i-channel of this theory; the Bogoliubov pp-ch&nnel 
is never relativistic. Consequently, it is justified to approximate the pairing 
force by the best currently available non-relativistic interaction: the pairing 
part of the Gogny-force. This is certainly more realistic than a one-meson 
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exchange force in the pairing channel, since this type of interactions have 
never been optimized for the description of pairing properties. 

2.3 Deformed Dirac-Hartree-Bogoliubov Equations 

The self-consistent solution of the Dirac-Hartree-Bogoliubov integro-differential 
eigenvalue equations and Klein-Gordon equations for the meson fields deter- 
mines the nuclear ground state. For systems with spherical symmetry, i.e. 
single closed-shell nuclei, the coupled system of equations has been solved 
using finite element methods in coordinate space [|, |^, [35], [36], and by 



expansion in a basis of spherical harmonic oscillator || [8], ^] . For deformed 
nuclei the present version of the model does not include solutions in coor- 
dinate space. The Dirac-Hartree-Bogoliubov equations and the equations 
for the meson fields are solved by expanding the nucleon spinors Uk(r) and 
Vfc(r), and the meson fields, in terms of the eigenfunctions of a deformed axi- 
ally symmetric oscillator potential ||28|| . For nuclei at the drip-lines, however, 



solutions in configurational representation might not provide an accurate 
description of properties that crucially depend on the spatial extension of 
nucleon densities, as for example nuclear radii. In less exotic nuclei on the 
neutron rich side, or for proton-rich nuclei, an expansion in a large oscillator 
basis should provide sufficiently accurate solutions. In particular, proton-rich 
nuclei are stabilized by the Coulomb barrier which tends to localize the pro- 
ton density in the nuclear interior and thus prevents the formation of objects 
with extreme spatial extension. 

The current version of the model describes axially symmetric deformed 
shapes. For this geometry the rotational symmetry is broken and the total 
angular momentum j is no longer a good quantum number. However, the 
densities and fields are still invariant with respect to a rotation around the 



symmetry axis, which is taken to be the z-axis. It is then natural to use 
cylindrical coordinates to represent the solutions of the RHB equations. In 
particular, the single-nucleon spinors U(r) and V(r) in (9) are characterized 
by the quantum numbers [Q, 7r,£ 3 ]. Q is the eigenvalue of the symmetry 
operator j z (the projection of the total angular momentum on the symmetry 
axis), 7r is the parity, and £3 the isospin projection. The two Dirac spinors 
U(r) and V(r) are defined as 



U(r,s,t) = 



and 



V(r,s,t) 



fu(r,s,t) 
igu(r,s,t) 



fv(r,s,t) 
i9v(r,s,t) 



fu(z,r ± ) e ^ n+1 ^ 
V ig^r^eW+m* J 



( fv(^r ± ) e ^-V2)^ ^ 
fy(z,r ± ) e ^ +l ' 2 ^ 
ig+^r^e^- 1 / 2 ^ 
V igv( z i r ±) e l(n+1/2)ip ) 



Xt 3 (t) (34) 



Xt 3 (t) (35) 



The basis expansion method starts from the deformed axially symmetric 
oscillator potential 

V osc (z,r ± ) = \muj 2 z z 2 + \Mu\r\. (36) 

Imposing volume conservation, the two oscillator frequencies hcu± and Hcu z 
can be expressed in terms of a deformation parameter (3 : 

hu z = huj exp(-W-^-/3 ) (37) 

V 47T 

1 / 5 

huj^ = huj exp(+-W— Po). (38) 

Z V 47T 

The corresponding oscillator length parameters are 



b z = \hr— and b± = (39) 



19 



Because of volume conservation, b 2 ± b z = 6jj. The basis is now specified by 
the two constants Tiujq and (3q. The eigenfunctions of the deformed harmonic 
oscillator potential are characterized by the quantum numbers 

| ol > — | n z ,n r ,mi,m s >, (40) 

where mi and m s are the components of the orbital angular momentum and 
the spin along the symmetry axis, respectively. The eigenvalue of j z , which 
is a conserved quantity in these calculations, is 

Q = mi + m s , (41) 

and the parity is given by 

7T = (-)»*+"*«. (42) 
The eigenvectors can be explicitly written as a product of functions 

* a (z,r ± ,<p,8,t) = <f> nz (z) <t>2(ri)^e im ^ X rns{s)Xt a {t) (43) 



with 



= $ Q (r, S ) Xta (0, (44) 



<PnM = %/WC)e-< 2/2 (45) 



N m i 

= f L ^f /J i:()l)e^ 1 (46) 

I 



and 



C = z/b„ V = r\jb\. (47) 

H n (() and L™{rj) denote the Hermite and associated Laguerre polynomials, 
respectively. The normalization constants are given by 



1 In I 
N n = . and N?' = J- r — -. (48) 
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For the nucleon spinors U(r) and V(r) in Eqs. fl3~3|) and (|35|), the following 
expansion in terms of the eigenfunctions is used 

fu(v)(r,s,t) = £ fu(y)a$ a (r,s)xtM ( 49 ) 

Ot-max 

9u(V)(r,s,t) = ^2 g u(v)a $ & (r,s)xti(t)- (50) 

a 

In order to avoid the occurrence of spurious states, the quantum numbers 
a max and a max are chosen in such a way, that the corresponding major quan- 
tum numbers N = n z + 2n p + mi are not larger than Np + 1 for the expansion 
of the small components, and not larger than Np for the expansion of the 
large components |f28fl . 

The solutions for the fields of massive mesons are also obtained with 
an expansion in a deformed oscillator basis. It is convenient to use the 
same deformation parameter (3q as for the fermion fields, and to choose the 
oscillator length bs = bo/\/2 : 

1 Nb 

<K*,r±) = -st^ 2 / 2 -"/ 2 E KnrK, fUC) y^LKn), (si) 

Og n z n r 

with ( = y2z/b z and 77 = 2r\jb\. Inserting this ansatz into the Klein- 
Gordon equation, an inhomogeneous set of linear equations is obtained 

N B 

E 7~tn z n r n z >n r > 4>n z 'n r ' = S n*n r ' (^) 
n z 'n/ 

with the matrix elements 

T~tn z n r n z 'n r ' = m t ^>n r n r l ^n z n z ' (53) 

+ ((2n 2 + 1) 5„ zn2 .' — n/^'nz+i — n 2 <5 nz „ z ' + i) <5 nrnr / 

+ "To" ((2«r + 1) <^n r n r ' + ^r'^n r 'ri r +l + «r ^n r n r '+l) ^n z n 2 ' 

Most of our calculations are performed by expansion in 12 oscillator shells 
for the fermion fields, and 20 shells for the boson fields. In order to test the 
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convergence, critical cases are also calculated with 14 shells for the fermion 
fields. 

2.4 Systems with Odd Number of Particles 

In the study of proton emitters nuclei with both even and odd number of 
protons (neutrons) must be calculated. The RHB wave function of a system 
with an odd number of particles is represented by the product state 
in the canonical basis. Therefore, for a nucleus with odd proton (neutron) 
number, one proton (neutron) level must be blocked. If both Z and N are odd, 
one proton level and one neutron level must be blocked. As described in Ref. 
[j3lH , a quasiparticle level k is blocked by exchanging the two eigensolutions 
(U k , V k ) and (V k , U k ) of the RHB-equations (|20]). The blocking procedure is 
performed in each step of the self-consistent iteration. The resulting density 



matrix and pairing tensor, in the oscillator basis, read |38| 

Paa' = Paa' + ^ akV a'k ~ ^ak^a'ki (54) 



K aa' — K aa' + ^akYa'k ^akJ^a'ki (55) 

where 

Paa' = (V*V T ) aal and n aa , = {V*U T ) aa , (56) 

are the corresponding matrices for the system without blocking. Each blocked 
level, of course, produces a different solution for the ground state. It might 
happen that several quasiparticle levels available to the odd proton (neutron) 
are found close in energy. This case is especially important for deformed nu- 
clei. In order to determine the self-consistent solution for the ground state, 
each available level has to be blocked in turn. The ground state then cor- 
responds to the self-consistent solution with lowest energy. Separation en- 
ergies and spectroscopic factors of proton emitters are calculated from the 
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self-consistent ground state RHB solutions for the parent system (odd Z) and 
daughter nucleus (even Z). Of course, both nuclei can have an odd number 
of neutrons, in which case calculations with blocked neutron levels must be 
performed. 

3 Proton drip-line and deformed ground-state 
proton emitters 

In the present work the relativistic Hartree-Bogoliubov theory is applied 
in the description of ground-state properties of deformed proton-rich odd 
Z nuclei in the region 53<Z<69. In particular, the RHB model is used 
to study the location of the proton drip-line, the ground-state quadrupole 
deformations and one-proton separation energies at and beyond the drip-line, 
the deformed single-particle orbitals occupied by the odd valence proton, and 
the corresponding spectroscopic factors. 

The input parameters of the RHB model are the coupling constants and 
the masses for the effective mean-field Lagrangian, and the effective interac- 
tion in the pairing channel. As we have done in most of our recent studies, 
we use the NL3 effective interaction |H| for the RMF Lagrangian. Properties 
calculated with NL3 indicate that this is probably the best effective interac- 
tion so far, both for nuclei at and away from the line of /^-stability. A very 
recent systematic theoretical study of ground-state properties of more than 
1300 even-even isotopes has shown excellent agreement with experimental 
data [fHJ. For the pairing field we employ the pairing part of the Gogny 
interaction with the parameter set D1S As we have already discussed 
in the previous section, this force has been very carefully adjusted to the 
pairing properties of finite nuclei all over the periodic table. 



23 



In Figs. |l| and |2] the one-proton separation energies 

S P (Z, N) = B(Z, N) - B(Z -1,N), (57) 

are displayed for the odd-Z nuclei 53 < Z < 69, as function of the number of 
neutrons. The model predicts the drip-line nuclei: 110 I, 115 Cs, 118 La, 124 Pr, 
129 Pm, 134 Eu, 139 Tb, 146 Ho, and 152 Tm. For nuclei with Z > 50, the su- 
perposition of the Coulomb and centrifugal potentials results in a relatively 
high potential energy barrier, through which the tunneling of the odd proton 
proceeds. For the ground-state proton emission to occur, the valence proton 
must penetrate the wide potential barrier, and this process competes with 
f3 + decay. The half-life of the decay strongly depends on the energy of the 
proton, i.e. on the Q p value, and on the angular momentum of the proton 
state. The Q p dependence has been nicely illustrated in Ref. [18|], where 



proton partial decay half-lives for 147 Tm have been calculated in the DWBA. 
It has been shown that, as the Q p value changes from 0.5 MeV to 2.5 MeV, 
the half-life is reduced by more than 22 orders of magnitude, from 10 10 sec 
to 10 -12 sec for the decay of the ld 3 / 2 state. This means that in many nuclei 
the ground-state proton emission will not be observed immediately after the 
drip-line, because for small values of Q p the total half-life will be completely 
dominated by j3 + decay. Only for higher values of Q p is proton radioactiv- 
ity expected to dominate over (3 + decay. By further decreasing the number 
of neutrons, the odd proton energy increases rapidly and this results in ex- 
tremely short proton-emission half-lives. For a typical rare-earth nucleus the 
Q p window in which ground-state proton decay can be directly observed is 



about 0.8 - 1.7 MeV fllSfl . An interesting point is that for nuclei at the 
drip-lines the precise values of the separation energies should also depend 
on the isovector properties of the spin-orbit interaction. It has been shown 
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that the relativistic mean-field model encloses a unique isovector dependence 
of the spin-orbit term of the effective single-nucleon potential 0, very dif- 
ferent from the one included in the more traditional non-relativistic models 
based on Skyrme ansatz. It would be particularly interesting to compare 
predictions of various models based on the framework of the self-consistent 
mean-field theory, as it was done, for example, in Ref. HI] for the structure 
of proton-drip nuclei around the doubly magic 48 Ni. 

In Figs. |l] and || we display the separation energies of nuclei beyond the 
drip-line. The energy window extends to include those nuclei for which a 
direct observation of ground-state proton emission is in principle possible 
on the basis of calculated separation energies. For the most probable pro- 
ton emitters, in Table we enclose the ground-state properties calculated 
in the RHB model. For each nucleus we include the one-proton separation 
energy S p , the quadrupole deformation f3 2 , the deformed single-particle or- 
bital occupied by the odd valence proton, and the corresponding theoretical 
spectroscopic factor. In order to determine the Nilsson quantum numbers 
and the spectroscopic factors, we analyze the wavefunctions of the parent 
nucleus with an odd proton number and that of the daughter nucleus with 
even proton number in the canonical basis. Because of the different mean 
fields and because of the blocking, the canonical basis-systems found in this 
way are not completely identical for the even and for the odd case. However, 
the wavefunctions and the occupation numbers v 2 are in most of the cases 
rather similar. The essential difference consists in the fact, that one of the 
pairs of the degenerate occupation numbers v 2 does not occur in the odd 
system. Instead of that we find the occupation numbers and 1. This level 
is obviously the blocked level fiQ in the canonical basis. The Nilsson quantum 
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numbers of the odd-proton level are determined by the dominant component 
in the expansion of this wavefunction in terms of the anisotropic oscillator 
basis (0). The spectroscopic factor of the deformed odd-proton orbital is 
defined as the probability that this state is found empty in the daughter 
nucleus with even number of protons. If one neglects the difference of the 
canonical basis for the even and the odd proton system, the spectroscopic 
factor of the level //o is calculated from the Eqs. (|30|) and ([3l]): 

= l(^ol«Iol^)| 2 = <• (58) 

where m w is determined in the canonical basis of the daughter nucleus with 
even proton number. 

The results of RHB calculations are compared with the predictions of the 
finite-range droplet (FRDM) mass model: the projection of the odd-proton 
angular momentum on the symmetry axis and the parity of the odd-proton 
state flp |]13|| , the one-proton separation energy fl3|| , and the ground-state 



quadrupole deformation [|T^]. The separation energies are also compared 
with available experimental data on proton transition energies. An excellent 
agreement is found between ground-state quadrupole deformations calculated 
in the two theoretical models. In particular, both models predict a transition 
from prolate to oblate shapes in the Tm isotopes. On the other hand, the 
predictions for the deformed orbitals occupied by the odd proton differ in 
many cases. For the I and Cs isotopes, the RHB calculations predict the 
3/2 + [422] and l/2 + [420] proton orbitals, respectively. Exactly the opposite 
ordering of orbitals is calculated in the macroscopic-microscopic model. The 
differences are also pronounced in 115,116 La, 130,131 Eu, and 145 Tm. In order 
to determine which is the last occupied proton orbital, the corresponding 
decay width for proton emission should be calculated and compared with the 
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experiment value. For spherical proton emitters this is a relatively simple 
task [|18| , but the extension of the formalism to deformed nuclei is rather 
involved. Calculations based on the coupled-channel formalism for deformed 
proton emitters have been reported for 109 I and 113 Cs [I5|, and more recently 
in Ref. [4~3| for 113 Cs. These calculations are not based on a fully microscopic 



and self-consistent description of proton unstable nuclei, but use a Woods- 
Saxon potential adjusted in such a way that the energy of the proton quasi- 
bound state coincides with the experimental transition energy. In particular, 
even the most realistic calculations do not include pairing correlations, or it 
is just the spectroscopic factor that is calculated in the simple BCS scheme. 



In Ref. [42] it was shown that the theoretical half-lives calculated using the 
proton orbital 1/2+ [420] and (3 2 « 0.05 - 0.10 for 109 I, and 3/2+[422] and 
(3 2 ~ 0.10 — 0.15 for 113 Cs, are consistent with the experimental values. 

The RHB and FRDM models predict similar values for the one-proton 
separation energies of proton emitters with Z > 59. Again, the differences 
are largest for I and Cs. The RHB predicts stronger binding for I, and 
weaker binding for Cs. When compared with available experimental proton 
energies, the FRDM model is very close to experiment for 109 I and 112 Cs, 
while the RHB model does better for 113 Cs. However, both models fail to 
reproduce the observed anomaly in the one-proton separation energies of 
112 Cs and 113 Cs. The experimental data indicate that 112 Cs is more bound, 
or better, less unbound. This discrepancy is easily explained. While in the 
odd-even nucleus 112 Cs the neutrons are all paired off, 113 Cs has an odd 
number of protons and an odd number of neutrons. Since N — Z is only 2, 
one expects a relatively strong interaction between the odd neutron and the 
odd proton. Compared to the odd-even system, the additional interaction 
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will increase the binding energy. Such an additional interaction could be 
represented, for instance, by a surface delta-force. Since the two mean-field 
models do not include any residual proton-neutron interaction, they cannot 
reproduce the inversion of the separation energies. The same effect could be 
also expected for the proton separation energies of 108 I and 109 I. Ground-state 
proton emission has not been observed in 108 I. The nucleus has been found to 
be an alpha emitter, with an upper limit for the proton decay branch of 1%. 
This limit implies a reduction in the Q p value of at least 220 keV compared 
to 109 I [JHJ, and the origin of this additional binding is again the residual 
interaction between the odd proton and the odd neutron. 

The theoretical separation energies are also compared with recently re- 
ported experimental data on proton radioactivity from 131 Eu, 141 Ho [H] . 
145 Tm ||44|| , 146 Tm [55], and 147 Tm JTjJ. The 131 Eu transition has an energy 
E p = 0.950(8) MeV and a half-life 26(6) ms, consistent with decay from ei- 
ther 3/2 + [411] or 5/2 + [413] Nilsson orbital. For 141 Ho the transition energy 
is E p = 1.169(8) MeV, and the half-life 4.2(4) ms is assigned to the decay of 
the 7/2~[523] orbital. The calculated RHB proton separation energy, both 
for 131 Eu and 141 Ho, is -0.9 MeV. In the RHB calculation for 131 Eu the odd 
proton occupies the 5/2 + [413] orbital, while the ground state of 141 Ho cor- 
responds to the 7/2~[523] proton orbital. This orbital is also occupied by 
the odd proton in the calculated ground states of 145 Tm, 146 Tm and 147 Tm. 
For the proton separation energies we obtain: —1.43 MeV in 145 Tm, —1.20 
MeV in 146 Tm, and -0.96 MeV in 147 Tm. These are compared with the 
experimental values for transition energies: E p = 1.728(10) MeV in 145 Tm, 
E p = 1.120(10) MeV in 146 Tm, and E p = 1.054(19) MeV in 147 Tm. When 



compared with spherical WKB or DWBA calculations [18], the experimental 
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half-lives for the three Tm isotopes are consistent with spectroscopic factors 
for decays from the hn/2 proton orbital. Though our predicted ground-state 
configuration 7/2~[523] indeed originates from the spherical hn/2 orbital, we 
find that these nuclei are deformed. 145 Tm has a prolate quadrupole de- 
formation 02 = 0.23; 146 Tm and 147 Tm are oblate in the ground-state with 
fa = —0.20 and fa = —0.19, respectively. Calculations also predict possible 
proton emitters 136 Tb and 135 Tb with separation energies —0.90 MeV and 
— 1.15 MeV, respectively. In both isotopes the predicted ground-state proton 
configuration is 3/2 + [411]. Another possible proton emitter is 130 Eu with 
separation energy —1.22 MeV and the last occupied proton orbital 5/2~[532] 
or 5/2 + [413]. In our calculation these two levels are almost degenerate, and 
the resulting binding and separation energies are identical. 

The calculated mass quadrupole deformation parameters for the odd-Z 
nuclei 53 < Z < 69 at and beyond the drip line are shown in Fig. |3|. While 
prolate deformations 0.15 < fa < 0.20 are calculated for most of the I and 
Cs nuclei, the proton-rich isotopes of La, Pr, Pm, Eu and Tb are strongly 
prolate deformed (fa ~ 0.30 — 0.35). By increasing the number of neutrons, 
Ho and Tm display a transition from prolate to oblate shapes. The absolute 
values of fa decrease as we approach the spherical solutions at N = 82. The 
quadrupole deformations calculated with the RHB /NL3 effective interaction 
are found in excellent agreement with the predictions of the macroscopic- 
microscopic mass model [12[ . 

A detailed analysis of single proton levels, including spectroscopic fac- 
tors, can be performed in the canonical basis which results from the fully 
microscopic and self-consistent RHB calculations. For the Eu isotopes this 
is illustrated in Fig. f|, where we display the proton single-particle energies 
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in the canonical basis as function of the neutron number. The dashed line 
denotes the position of the Fermi level. The proton energies are the diagonal 
matrix elements of the Dirac Hamiltonian hp in the canonical basis. The 
phase-space of positive-energy states should not be confused with the con- 
tinuum of scattering states which asymptotically behave as plane waves. The 
RHB ground-state wave function can be written either in the quasiparticle 
basis as a product of independent quasi-particle states, or in the canonical 
basis as a highly correlated BCS-state. In the canonical basis nucleons oc- 
cupy single-particle states. The canonical states are eigenstates of the RHB 
density matrix. The eigenvalues are the corresponding occupation numbers. 
In particular, we notice that for the proton emitter 131 Eu, the ground-state 
corresponds to the odd valence proton in the 5/2 + [413] orbital. For 130 Eu the 
states 5/2 + [413] and 5/2 _ [532] are almost degenerate, and this is an example 
of a situation in which only the comparison of calculated and measured half- 
lives can decide which is the last occupied proton orbital. Both for 130 Eu and 
131 Eu the macroscopic-microscopic model predicts the last proton orbital to 
be fi 71 " = 3/2 + , which corresponds to the Nilsson orbital 3/2 + [411]. In the 
energy diagram of Fig. f|the states 5/2 + [413] and 5/2 _ [532] are closer to the 
Fermi level for all calculated Eu isotopes. 

In Fig. [5] we plot the proton single-particle energies in the canonical basis 
for the nuclei at the drip-line: 118 La, 124 Pr, 129 Pm, 134 Eu, 139 Tb. The levels 
are shown as functions of the number of protons, and again the dashed line 
denotes the position of the Fermi level. From this diagram one can easily 
deduce which levels are most likely to be occupied by the odd-proton, at 
and beyond the drip-line. Figure should be compared with Table [I] for the 
odd-proton orbitals occupied in the proton emitters. Of course, by further 
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decreasing the number of neutrons, i.e. by going beyond the drip-line, the 
proton levels and the Fermi level are shifted upwards in energy. However, as 
can be seen from Fig. ^ for the Eu isotopes, this increase in energy is very 
smooth and dramatic changes should not be expected for the occupation of 
the proton orbitals beyond the drip-line. 

In many applications of the HFB and RHB models it has been emphasized 
how important is the self-consistent treatment of pairing correlations in drip- 
line nuclei. In order to illustrate the relative contribution of the pairing 
correlations to the binding energy, for the drip-line nuclei 51 < Z < 69, in 
Fig. | we display the average values of the proton pairing gaps for occupied 
canonical states 

< A p >= Eln ' |A, TK (59) 

where vf Qn i are the occupation probabilities of the proton orbitals, and A^j 
are the diagonal matrix elements of the pairing part of the RHB single- 
nucleon Hamiltonian in the canonical basis. The plotted values of < A p > 
correspond to the RHB self-consistent solutions for the ground-states of the 
drip-line nuclei. The pairing correlations vanish for Sb and I. The average 
proton pairing gap is 1 — 1.5 MeV for drip-line nuclei between 115 Cs and 
139 Tb, and it increases to more than 2 MeV for 146 Ho and 152 Tm. 



4 Conclusions 

The relativistic Hartree-Bogoliubov theory has been used to study the ground- 
state properties of deformed proton emitters with 53< Z < 69. The principal 
advantage of working in the relativistic framework is its ability to describe 
all the nuclear mean fields in a unified and fully self-consistent way. In 
particular, the spin-orbit term of the effective single-nucleon potential is au- 
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tomatically included in the relativistic mean-field model. In contrast to non- 
relativistic models, the strength and isospin dependence of the spin-orbit in- 
teraction do not require the introduction of additional free parameters. With 
the addition of the pairing force, the relativistic Hartree-Bogoliubov theory 
provides a unified and self-consistent description of mean-field and pairing 
correlations. This theoretical framework is especially important for applica- 
tions in exotic nuclei far from the valley of /3-stability. The model parameters 
are the coupling constants and meson masses of the mean-field Lagrangian, 
and the interaction in the pp channel. The NL3 parameter set has been used 
for the mean-field Lagrangian, and pairing correlations have been described 
by the pairing part of the finite range Gogny interaction D1S. This particular 
combination of effective forces in the ph and pp channels has been used in 
most of our applications of the RHB theory. The calculated ground-state 
properties have been found in excellent agreement with experimental data, 
not only for spherical /3-stable nuclei, but also for nuclear systems with large 
isospin values on both sides of the valley of /3-stability. No attempt has been 
made to adjust the model parameters to the specific properties of nuclei 
studied in this work. 

The model has been used to study the location of the proton drip-line, the 
ground-state quadrupole deformations and one-proton separation energies at 
and beyond the drip line, the deformed single particle orbitals occupied by the 
odd valence proton, and the corresponding spectroscopic factors. The results 
of fully self-consistent calculations have been compared with predictions of 
the finite range droplet mass model, and with very recent experimental data 
on ground-state proton emitters. 

The calculated ground-state mass deformation parameters have been found 
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in excellent agreement with values predicted by the macroscopic-microscopic 
FRDM model calculations. In particular, both models predict a transition 
from prolate to oblate shapes in the Tm isotopes. In many cases, however, 
different single particle orbitals occupied by the odd proton are predicted by 
the two models. Correspondingly, in these cases relatively large differences for 
the one-proton separation energies have been calculated. These divergences 
might be related to the different descriptions of the effective spin-orbit po- 
tential. In comparison with recent experimental data on ground-state proton 
emitters, the calculated one-proton separation energies are in excellent agree- 
ment with transitions energies in 113 Cs, 131 Eu, 141 Ho, 146 Tm, and 147 Tm. For 
145 Tm, however, the calculated and experimental proton energies differ by 
almost 300 keV. Both the RHB and FRDM models fail to reproduce the 
observed anomaly in the proton energies of the ground state emitters 112 Cs 
and 113 Cs. The anomaly can probably be related to the short range proton- 
neutron interaction, not included in the mean-field models. Summarizing, we 
have found a remarkable agreement between results of RHB calculations and 
very recent experimental data on deformed ground-state proton emitters. It 
should be emphasized that these results have been obtained with the stan- 
dard RHB parameterization, i.e. no new parameters have been introduced 
to specifically describe the proton emitters. In addition, the spectroscopic 
factors which correspond to the fully self-consistent RHB solutions will be 
useful in calculations of proton partial decay half-lives. 
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Figure 1: Predictions of the RHB model (NL3 + D1S) for the one-proton 
separation energies of odd-Z nuclei 61 < Z < 69, at and beyond the drip-line. 



Figure 2: Same as in Fig. [T], but for odd-Z isotopes 53 < Z < 59. 



Figure 3: Self-consistent ground-state quadrupole deformations for odd-Z 
nuclei 53 < Z < 69, at and beyond the proton drip-line. 



Figure 4: The proton single-particle levels for the Eu isotopes. The dashed 
line denotes the position of the Fermi level. The energies in the canonical 
basis correspond to ground-state solutions calculated with the RHB model 
(NL3 + D1S). 



Figure 5: The proton single-particle levels for the isotopes at the proton 
drip-line 57 < Z < 65. The energies in the canonical basis are plotted as 
functions of the proton number Z, and the dashed line denotes the position 
of the Fermi level. 



Figure 6: Average proton pairing gaps < A p > for the isotopes at the proton 
drip-line: U0 I, 115 Cs, 118 La, 124 Pr, 129 Pm, 134 Eu, 139 Tb, 146 Ho, and 152 Tm. 
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Table 1: Odd-Z ground-state proton emitters in the region of nuclei with 
53 < Z < 69. Results of the RHB calculation for the one-proton Separation 
energies S p , quadrupole deformations /3 2 , and the deformed single-particle 
orbitals occupied by the odd valence proton, are compared with predictions 
of the macroscopic-microscopic mass model, and with the experimental tran- 
sition energies. All energies are in units of MeV; the RHB spectroscopic 
factors are displayed in the sixth column. 
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